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Abstract. Recent X-ray observations reveal an increasing number of X-ray sources in nearby galaxies exceeding 
luminosities of Lx ^ 2 • 10®®erg s“^. Assuming isotropic e mission, the Eddington limit sugg ests a population of 
intermediate-mass black holes of M, ';$> IOMq. However, Vlarkoff, Falcke, & Fender (2001) proposed that jets 
may be contributing to the X-ray emission from X-ray binaries (XRBs), implying that some X-ray sources may 
be relativistically beamed. This could reduce the required black hole masses to standard values. To test this 
hypothesis, we investigate a simple X-ray population synthesis model for X-ray point sources in galaxies with 
relativistic beaming and compare it with an isotropic emission model. The model is used to explain a combined 
data set of X-ray point sources in nearby galaxies. We show that the current distributions are consistent with 
black hole masses Mq ^ 10 and bulk Lorentz factors for jets in microquasars of yj ~ 5. Alternatively, intermediate 


mass black holes up to 1000 Mq are required which are distributed in a powerlaw with roughly 


M" 
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1. Introduction 

In recent years X-ray observations have revealed several 
off-nucleus ultra-luminous X-ray sources (ULXs) in the lu- 
minosity range 10^^ — 10^°erg s~^ within nearby galaxies 


(La Parola et al. 200 it 

Mizuno, Kubota, & Makishima 

2001; 

Bauer et al. 2001 

; IColbert & Mushotzky 19991). 


this is difficult to achieve with pu re disk models. Recently 
Markoff, Falcke, fc Fender (2001 ) suggested that the spec¬ 
trum of some XRBs could be explained by a coupled 
disk/jet model, where some of the X-ray emission is pro¬ 
duced by synchrotron and inverse-Compton radiation in 
a jet. This emission would naturally be relativistically 
beamed. Mirabel fc Rodriguez (199^ ) (see also Reynold; 


The E jldington limit for an accreting object with mass out that a number of nearby 

% 


M is LEdd ~ 1.25 • 10“®j^6fg S which iirlplles that 


these sources are super-Eddington for stellar mass ob¬ 
jects. Some ULXs show spectral transitions from a soft 
spectrum to a hard power law and luminosity variability 
(e.g. Mizuno, Kubota, fc Makishima 2001 ; Kubota et al. 
2001 ), ruling out supernova remnants and supporting the 
idea that ULXs can be attributed to accreting black holes. 
To achieve the observed X-ray luminosities with isotropi¬ 
cally radiating accretion disks should require a population 
of intermediate-mass black holes of 50 — SOOMq. As dis¬ 
cussed in Kubota et al. (2001), however, the measured 
inner-disk temperatures of ULXs (Tin = 1.0 — 1.8 keV) 
are too high for these masses and there is no established 
formation scenario for such intermediate-mass black holes. 

The problems with i sotropic emission models have al¬ 
ready been discussed by King et al. (2001 ), where the au¬ 
thors propose some form of anisotropic emission as an al¬ 
ternative. A beaming factor of ten already reduces the 
required mass of the black holes to expected values, but 
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galaxies should host microblazars - microquasars with rel¬ 
ativistically beamed jets pointed towards the observer. We 
will here investigate whether such populations of microb¬ 
lazars or intermediate mass black holes can indeed explain 
current data on ULXs and constrain the basic parameters 
required for these models. 

2. Simple models 

2.1. The Jet/Disk Model 

Black hole candidate XRBs can exist in several states, the 
two most distinct of which are a high/soft state where 
the observed spectrum is soft and thermally-dominated, 
and a low/hard state dominated by a non-thermal hard 
power law spectrum (e.g. Nowak 1995| ). Which state an 
XRB is in seems to depend on the accretion rate. One sce¬ 
nario for the evolution of XRBs is that the inner part of 
the accretion disk consists of an optically thin, advection- 
dominated accretion flow (ADAF) existing up to a tran¬ 
sition radius where the accretion flow turns into a stan- 
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dard ( Shakura fc Sunyaev 1973 ) optically thick disk ( Esin 
McCliitock, & Narayan 1997). The low/hard state seems 
to be accompanied by persistent radio jets with optically 
thick synchrotron emission extending up to the optical 
and near-infrared ( Fender 20011) . Synchrotron and inverse 
Compton emission from the jet could also be produced 


Hence, we use the following simple parameterization 
for the soft X-ray luminosity of accretion disk and jet: 


in the X-rays for low/hard and high/soft states (Markoff, 
Falcke, & Fender 2001; Markoff et al. 2001). 

For a self-consistent population synthesis model, we 
will have to take the jet and the disk separately into ac¬ 
count. For simplicity we make some assumptions in order 
to calculate the relative importance of the two processes 
for the overall distribution of XRBs. These are: 

— We only consider two main populations: neutron stars 
of mass 1 . 4 M 0 and black holes within a mass range 
of 5-15 Mq, where we fix the amount of active black 
holes ( > 5- lO^^erg s~^) to 13% of the number of active 


Tdisk — 


Tjet — 




if M<Mc 
eMc^ if Me < M < iHEdd 


if M<Mc 
rjeMcP if Me < M < MEdd 




( 1 ) 


neutron stars (e.g. Tanaka fc Lewin 1995 ). 

— While the spectral states are defined only for black 
holes, for simplicity we apply them also to neutron 
stars, and consider only the low/hard and soft/high 
states. 

— The probability that a given XRB has the accretion 
rate M is given by W{M) ~ M^ which we assume as 
a power law with a cutoff representing the Eddington 
limit. 

— The distribution of accretion rates and basic jet pa¬ 
rameters are assumed to be identical for neutrons stars 
and black holes. 

— Soft X-ray emission is produced by an isotropically ra¬ 
diating disk and a relativistically beamed jet as dis¬ 
cussed below. 

The crucial point for a simple population synthesis 
model including jet emission is how a specific accretion 
rate translates into an X-ray luminosity. Here we assume 
that the transition between the low and the high state 
happens at a critical accretion rate Me and the lumi¬ 
nosities scale as follows: above Me the disk luminosity 
increa ses linearly with M as expected for a standard ac- 


In the following we set the radiative efficiency of the stan¬ 
dard accretion disk to the canonical value of e = 0.1. For a 
given mass M the parameter MEdd has been chosen such 
that the luminosity of the disk and the jet integrated over 
all angles is equal to Tfidd- For simplicity, the mass distri¬ 
bution of black holes is given by dN/dM = V{M) =const. 

For a bulk Lorentz factor of the jet of yj > 1 the jet 
emission depends on the angle to the line of sight as given 
by Lind fc Blandford (1985| ). If the emission in the rest 
frame of the jet follows a power-law with spectral index 
a, the observed emission is proportional to where 

the Doppler factor S = ^.(^i_pcose) ■ probability of 
seeing an object with an emission exceeding L when in 
the rest frame the jet emits Lioc is: 


P(L, Lioc) = 



- 1 


where Tmax = (5^'''“(0 = 0 )Lioc is the maximum possible 
emission. To derive this we only consider the jet compo¬ 
nent pointing towards us and then integrate over all incli¬ 
nation angles. Since we only discuss jets with yj > 2, the 
emission of the counter-jet is largely negligible. 

With this parameterization, the contribution of a sin¬ 
gle population of XRBs at a given accretion rate and mass 
has three parameters (yj, M^o, 77 ), which are reasonably 
well constrained by the underlying models. Observations 
of microquasars show that the typical Lorentz factors are 


cretioni disk. Below Me. the disk luminosity increases with 


in the range yj ~ 2 —5 (Mirabel & Rodriguez 1999; Fender 


M^ as expected for optically thin ADAFs ( Narayan fc 
Yi 199|; for a constant a-parameter). Assuming that the 
jet power scales linearly with M, the optically thick jet 
synchrotron emission will scale roughly as Tx.jet o<^ 
(Falcke & Biermann 1995; Falcke fc Biermann 1999| ). For 
simplicity we assume that this scaling also holds for the 
optically thin flux. At high accretion rates this scaling 
must break down when a significant fraction of the jet 
power is radiated away. In this phase the radiated power 
can only increase linearly with jet power. In the jet model 
of Markoff, Falcke, fc Fender (2001D and Markoff et al. 
( 2001 |)Jthis happens in the high state, roughly at M > Me, 
where the jet is inverse-Compton cooled (radiating soft X- 
rays) by the accretion disk. However, since models for the 
contribution of jets to the high state of XRBs are not yet 
very well developed, we simply fix the luminosity of the 
jet at Ljet = ??Tdisk at Me, where 77 is a free parameter. 


et al. 1999 ). The typical values for the critical accretion 
r ate discussed in the literature are around Me ^ 0.1 
( Narayan fc Yi 1995 ) and we keep this parameter fixed. 
The jet efficie ncy can in principle be fai rly high, but prob¬ 
ably r] ^ 0.3 ( Falcke fc Biermann 1995 & 1999). 

At high accretion rates, the jets emits a factor 77 less 
radiation in its rest frame compared to the disk, but beam¬ 
ing will lead to an amplification of the jet with respect to 
the disk for small inclination angles. Beaming is strongest 
for sources whose jets point within the beaming cone with 
half-opening angle of l/yj towards the observer. For ex¬ 
ample, already yj = 3 will beam the jet component in a 
fraction of 5.5% of the binaries by a factor of 20. This is 
more than enough to make up for the less efficient emis¬ 
sion mechanism. For yj = 5 a fraction of 2% binaries are 
beamed by a factor of 77. Therefore, in the low-luminosity 
regime (<< cqMe) jets should dominate (because of low 
radiative efficiency of ADAFs), and again dominate in the 
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super-Eddington regime due to beaming. In the interme¬ 
diate regime up to LEdd, the disk will be more prominent. 

Now we can calculate a synthetic log N - log L distri¬ 
bution for our model. The emission from the disk and the 
jet is described by Eq. To reach a given luminosity L, 
only the difference L — Ldisk has to be reached by the jet 
due to boosting. The estimated number of sources with a 
luminosity greater than L is given by: 

N{L) = ^ M /dM / dMVi{M)Wi{M) ■ 

i=N,B •' •' 

-P(T-Tdisk(M),Ljet(M)) (2) 

where the sum goes over the two populations. 


2.2. The Disk-Only Model 

Of course, with the model discussed above we can also 
investigate the alternative scenario, that the soft X-ray 
emission originates only from the accretion disk, by setting 
77 = 0. In this case we have to leave the upper end of the 
black hole mass distribution and its power law index as a 
free parameter to obtain the high luminosities observed, 
and we have dN/dM = V{M) cx M‘^. 

With the isotropic disk emission the estimated number 
of sources with a luminosity greater than L is given by: 

N{L) = ^ M / dM / dMV*(M)W*(M) • 

_7\r D 


i^N,B 

0(L-Ldisk(M)) 

where 0 is the usual step-function. 


( 3 ) 


3. Data 

To put meaningful constraints on a beaming model it is 
essential to compare the low-luminosity (un-beamed) par¬ 
ent population with the high-luminosity (beamed) popula¬ 
tion. However, the published X-ray population of a single 
galaxy has only marginal statistics in the high-luminosity 
regime. To get a more general X-ray population in the low- 
luminosity regime up to « 5 • 10^®erg/s, we combine data 
from the galaxies MIOI , M3I and M82 ([Pence et al. 2001 


Di Stefano et al. 2001; [Griffiths et al. 200C| ). These are 
three close [D < 10 Mpc) galaxies with good published 
Chandra data. To get better statistics at higher luminosi- 
ties, we used the l uminosity function compiled by Roberts 
fc Warwick (200C| ) from ROSAT data of 49 spiral galaxies 
from the XHFS-sample. The host galaxy types are spiral 
with the exception of the irregular M82. However this does 
not seem to make a significant difference (see Fig. 1 & 2). 

To avoid incompleteness near the detection limit, for 
each data set we only use X-ray sources with a luminosity 
of ten times the respective detection threshold. Chandra 
has a different bandpass (0.3-10 keV) than ROSAT (0.1- 
2.2 keV), so we extrapolated the ROSAT-luminosities to 
the 2-8 keV band. If the photon index is not fitted directly 
we used a common powerlaw with Tph = 1.7. In some cases 
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Fig. 1. Comparison of our model of the luminosity func¬ 
tion with the data. The parameters are yj = 5 , ry = 0.3, 
Me = 0.1, ^ = 1.4. Also shown are the individual contri¬ 
butions of the disk and the jet for neutron stars and black 
holes. 


different values for Nh were used which we did not correct. 
According to Di Stefano et al. ( Di Stefano et al. 2001 ), 
variations of Nh — 0.6 — 1.5 • 10^^/cm^ and Tph = 1.2 — 2 
give differences in luminosity of about 20 % which are not 
really significant in the log(N)-log(S) plots shown here and 
should be statistically distributed (for ROSAT data the 
differences are higher). 

As a reference galaxy we take MlOl, to which we 
scale the populations of the other galaxies (i.e., the to¬ 
tal number of sources in the overlapping luminosity bins). 
Combining these data sets assumes that the overall shape 
of the luminosity distribution is roughly universal. Clearly, 
the overall number of XRBs in each galaxy can depend 
strongly on the age of recent star formation, but the aver¬ 
age slope of the luminosity function should be less sensi¬ 
tive to this. Since M82 is irregular and has a much higher 
star forming rate than MlOl or M31, we also show the 
data excluding M82, which is not significantly different. 

To calculate errors we assume a standard deviation 
from the ’general population’ of ^/N where N is the num¬ 
ber of detected sources, and use normal error propagation. 
Because we are showing a cumulative distribution, the er¬ 
rors for each point are not independent. 


4. Results 

To compare our simple model with the data we evaluate 
the integrals in Eqs. (^&^) numerically. The absolute nor¬ 
malization and the parameter which are entirely free, 
have been fit to the data at Lx < lO^^erg s“^. We obtain a 
best-fit value of the accretion rate index ^ = 1.4 (note that 
the luminosity scales as in this regime). As we only 
model neutron stars and black holes, the model fits could 
be affected at lower luminosities by other source popula¬ 
tions like accreting white dwarfs and supernova remnants. 
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Fig. 2. Model with intermediate-mass black holes up to 
IOOO-M 0 and a mass distribution with power law index 
C = 2. 


Fig. 0 shows the result for our best-fit jet/disk model 
which requires 7 j = 5 and 77 = 0.3, for the combined data 
set discussed above. The Eddington limit for black holes 
(limited to M < 15Mq) and neutron stars shows up as 
breaks at the respective luminosities, and jets with 7 j = 5 
are able to produce emission up to 10 ^°erg s“^ in signif¬ 
icant numbers. The model is most sensitive to 7 j and rj. 
Because the high luminosity domain depends linearly on rj 
while its dependence on 7 j goes as 7 ?'^, a slight decrease of 
7 j can be compensated by an increase of rj and vice versa. 
For 7 j = 5.8 or 7 j = 7.5 we can find rj = 0.2 or 77 = 0.1, but 
the fit gets progressively worse at higher Lorentz factors. 
Demanding 77 ^ 0.3 for the radiative efficiency of the jet 
sets a lower limit for 7 j ^ 5. 

For the disk-only model, the sensitive parameters are 
the power law indices of the accretion rate and the mass 
distribution of the black holes. To fit the data, a mass 
distribution index of ~ 2 is needed. The index of the ac¬ 
cretion rate is the same as before (1.4), because the lower 
luminosities in the jet/disk model are also dominated by 
the disk. To explain the most luminous sources, the up¬ 
per end of the black hole mass distribution must be ex¬ 
tended at least up to 1000 Mq. The results of the fit are 
shown in Fig. 0. For the lower luminosities the neutron 
stars dominate the luminosity function, while the black 
holes dominate at higher luminosities. 


5. Summary and Discussion 

Using two very simple models for the evolution of XRBs, 
we calculate the luminosity distribution of X-ray point 
sources i n nearby galaxies. We cons ider a jet/disk model 
based on Falcke fc Biermann (199^ ) and Markoff, Falcke, 
fc Fender (2001 ), which can give rise to relativistically 
beamed emission from microblazars. Alternatively we also 
consider a purely isotropically radiating disk model. 

Both models can in principle reproduce a combined lu¬ 
minosity function compiled from X-ray point source cat¬ 


alogs of three close galaxies and the XHFS spiral galaxy 
sample. However, as expected, the isotropic disk model re¬ 
quires a mass distribution of black holes extending out to 
IOOOM 0 to explain the ULXs. On the other hand, a rel¬ 
ativistic jet/disk model can fit the data with stellar mass 
black holes, if X-ray emitting jets with Lorentz factors 
7 j ~ 5 are present in XRBs. In addition, a fraction of 
77 = 10 — 30% of the total soft X-ray emission has to 
come from the jet rather than the accretion disk for an 
un-beamed XRB in the high state. This requires rather 
powerful jets but is not completely unreasonable. If only 
a fraction of the XRBs have relativistic jets, a slightly 
higher Lorentz factor or jet efficiency is needed. Boosting 
a 10 mJy Galactic XRB by a factor ^ 10^ (for 7 ^ 5 ) and 
placing it at D ~ 3Mpc would yield only a faint 10 nJy 
source and make radio detections difficult. 

With the current statistics it is not possible to dis¬ 
tinguish between the two different models, but it seems 
that microblazars provide at least a sensible alternative to 
the often discussed intermediate mass black hole scenario. 
Monitoring the spectral variability of the most luminous 
sources and further developing the XRB jet model should 
eventually help to disentangle the two scenarios. 
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